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o !■ The ~ 10 4 y old Vela Pulsar represents the bridge between the young Crab-like 

q . and the middle-aged rotation powered pulsars. Its multiwavelength behaviour is 

+3 , due to the superposition of different spectral components. We take advantage of 

the unprecedented harvest of photons collected by XMM-Newton to assess the 
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Vela Pulsar spectral shape and to study the pulsar spectrum as a function of its 
rotational phase. In order to fully exploit the data collected by XMM-Newton 



on the Vela pulsar, we had first to discriminate the pulsar emission from that of 
the bright surrounding nebula. To this aim, we used the Chandra/HRC surface 
brightness map of the nebula, coupled with the most accurate calibration of the 
EPIC point spread function. This procedure made it possible to assess the pulsar 



spectral shape, disentangling its thermal component from the non thermal one. 
As for the middle-aged pulsars Geminga, PSR B0656+14 and PSR B1055-52 (the 
"Three Musketeers"), the phase-integrated spectrum of Vela is well described by a 
three-component model, consisting of two blackbodies (Tbb = 1.06 ±0.03 x 10 6 K, 
R hh = 5.1+0-4 km? Tbb = 2.16±g;gf x 10 6 K, R BB = 0.73±g;g| km) plus a power-law 
(7 = 2.2lo 3)- The relative contributions of the three components are seen to vary 
as a function of the pulsar rotational phase. The two blackbodies have a shallow 
~ 7 — 9% modulation. The cooler blackbody, possibly related to the bulk of 
the neutron star surface, has a complex modulation, with two peaks per period, 
separated by ~ 0.35 in phase, the radio pulse occurring exactly in between. 
The hotter blackbody, possibly originating from a hot polar region, has a nearly 
sinusoidal modulation, with a single, broad maximum aligned with the second 
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peak of the cooler blackbody, trailing the radio pulse by ~ 0.15 in phase. The non 
thermal component, magnetospheric in origin, is present only during 20% of the 
pulsar phase and appears to be opposite to the radio pulse. XMM- Newton phase- 
resolved spectroscopy unveils the link between the thermally emitting surface of 
the neutron star and its charge-filled magnetosphere, probing emission geometry 
as a function of the pulsar rotation. This is a fundamental piece of information 
for future 3-dimensional modeling of the pulsar magnetosphere. 

Subject headings: pulsars:general - pulsars:individual (Vela) - stars: neutron - 
X-rays: stars 



1. Introduction 



The Vela Pulsar is one of the best scrutinized neutron stars. It was the first pulsating 
radio source observed in the southern hemisphere (ILarge et al.lll968l ) and the swing of the 
polarisation ve ctor during the radio pulse provided evidence for a rotational origin of the 
radio emission (jRadhakrishnan et al.lll969l ). Moreover the position of the pulsar, close to the 
centre of the Vela Supernova Remnant (hereafter SNR) confirmed the association of pulsars 
with rotating neutron stars born from massive stars' collapse. 



About ten years after the radio discovery, IWallace et al.l (119771 ) detected a puls ating op- 
tical source of V ~ 23.6: at least four peaks are present in the optical (IGouiffes 19981) and UV 



(Ro mani et al.l 120051 ) light curves, as well as in the hard X-ray energy range (IHarding et al. 



20021 ). HST observations (jCaraveo et al.l 120011 ) of the optical counterpart have allowed the 
direct measurement of the pul sar distance (294^^ pc), later confirmed and refined through 
radio observations (287^ pc iDodson et al.l 120031). PSR B0833-4 5 is also a bright 7-ray 



source: pulsations were detected by SAS-2 ([Thompson et al. 



975), COS-B (Bennett et al. 



19771 : iKanbach et all Il980l : iGrenier et al.l Il988f ) and CGRO. iKanbach et all (J199J), using 
EGRET, were able to perform spectroscopic analysis of different phase intervals of the dou- 
ble peaked 7-ray pulsar emission. 



ROSAT detection of soft X-ray pulsations (jOegelman et al.lll9931 ) was the last piece of 
the Vela multiwavelength puzzle. It turned out to be a difficult observation since X-rays 
from the neutron star are embedded in the bright Pulsar Wind Nebula (hereafter PWN), 



1 Based on observations with XMM-Newton, an ESA science mission with instruments and contributions 
directly funded by ESA member states and the USA (NASA). 
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located near the centre of the Vela SNR. A blackbody with temperature ~ 1.5 — 1.6 x 10 6 K 
and a radius of 3 — 4 km could describe ROSAT spectrum while the nebular emission, clearly 
non-thermal, could be ascribed to synchrotron emission originated from the interaction of 
the high-energy particle pulsar wind with the interstellar medium. Chandra observations 
of the Vela pulsar provided high resolution images of the X-ray nebula surrounding the 
neutron star. The nebula tu r ns out to be quite complex, with spectacular arc and jet-like 
features (IHelfand et al.l 1200 ll ; iPavlov et al.l 120031 ) reminiscent of those observed around the 
Crab. The thermal nature of the pulsar spectrum inferred from ROSAT data was confirmed. 
While high resolution spectra f ailed to show abso rption features, a non-thermal harder tail 
was found i n AC IS-S spectra (IPavlov et al.l 120011 ) . Analysing XMM-Newton observations 



Mori et al.l (120041 ) found the same thermal radiation. However, in their preliminary analysis, 
they studied only the pulsar phase-averaged thermal emission below ~ 1 keV. Since thermal 
emission is a very d istinctive character of mid dle aged pulsars (Geminga, PSR B0656+14 
and PSR B1055-52, feecker fc Truemperlfl997L Vela's (r « 11 400 y ; E « 7 x 10 36 erg s" 1 ) 
spectral properties make it more similar to older specimen than to younger ones. In this 
paper, we will provide further evidence in favor of this classification, taking advantage of 
XMM-Newton statistics to perform phase-resolved spectroscopy of the Vela pulsar. 



2. The Data 

We will use data collected by XMM-Newton and Chandra, exploiting the characteristics 
of both observatories. While XMM-Newton observations provide an unprecedented harvest of 
photons from the Vela pulsar and its PWN, the instrument point spread function (hereafter 
PSF, 6.6" FWHM), however, does not allow to disentangle the pulsar emission from the 
extended one. Accounting for the bright PWN contribution is mandatory to unveil the PSR 
emission properties. Thus, the high resolution image obtained by the Chandra observatory 
is crucial to have a clear view of the Vela pulsar and its sorrounding nebular emission. 



2.1. XMM-Newton EPIC Data 

The X-ray Multi Mirror Mission observed PSR B0833-45 on 2-3 December 2000. Owing 
to the source high flux, to reduce photon pile-up pn and MOS2 were used in "Small Window" 
mode, while MOS1 operated in "Large Window" . In the XMM-Newton Science Archive there 
are two data sets (id 0111080101 and 0111080201), separated by about 1.5 h, that cover a 
time span of (respectively) 41 100 s and 61 800 s. The raw data were first processed with 
standard pipeline tasks of the XMM Science Analysis Software (SASv6.5.0) and then the soft 
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proton flares were removed using the standard prescription of iDe Luca fc Molendil (l2004j ). 
The two data sets were merged using the SAS task merge. Table [I] summarise Good Time 
Intervals and Total Exposure time for the three EPIC detectors. 

Fig. [1] shows an image of the entire field of view, extracted using pn data: the compact 
nebula surrounding the Vela pulsar, with its peculiar shape, clearly stands out. The XMM- 
Newton telescopes angular resolution, however, does not allow to resolve the PSR from the 
PWN. Indeed no spatial selection can disentangle the point source from the diffuse emission. 



2.2. Chandra HRC and ACIS Data 



To get a clear view of the spatial distribution of the bright nebular emission, we analysed 
the two Chandra HRC-Imaging archive observations of the Vela pulsar, performed after the 
January 2000 glitch. A third observation, taken on January 2002 was also available. All 
these observations were reprocessed with CI AO0 v3. 2, applying aspect solution § and degap 
correction^, as well as reducing tap-ringing distortion^. The resulting image (Fig. [2]) shows 
the well known structure of the Ve la PWN: the two "arcs" and the "jet /co unter-jet" feature 
protruding from the point source (jlielfand et al. 2001 : Pavlov et al. 2003 ). 



3. Spectral Analysis 

To perform the spectral analysis of the XMM-Newton data on the Vela pulsar, we exploit 
both the spectral and temporal resolution of the EPIC/pn camera. To disentangle the PSR 
photons from the PWN ones, we take advantage of their different space distribution: while 
the PSR photons follow the instrument PSF, the nebular ones do not. 

To estimate the nebular contribution in the EPIC/pn event file, 

• we extract spectra from concentric annular regions of increasing radii; 

• we determine the PSR contribution in each spectrum as a function of the PSF and of 
the Encircled Energy Fraction (EEF); 

2 Chandra Interactive Analysis of Observations (CIAO), http://cxc.harvard.edu/ciao/ 
3 http : / / cxc . harvard . edu/ ciao/threads/ arcsec_correct ion/ 
http : / / cxc . harvard . edu/ciao/threads/hrci_degap/ 
D http : / / cxc . harvard . edu/ciao/threads/hrc^ampsf / 
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• we determine the PWN contribution in each spectrum using a surface brightness map 
of the nebular emission derived from Chandra data; 

• we fit all the spectra at once with a two component (PSR + PWN) model whose 
normalisation ratios are fixed following the two previous steps. 

In the following section we will describe in detail each step. 



3.1. EPIC/pn Phase- Averaged Spectra 

We extract spectra from the merged EPIC/pn event file, with FLAG and PATTERN 
0-4, from a 10" radius circle centered on the pulsar position as well as from three annular 
(10"-20", 20"-30", 30"-40") regions (see Fig. [p. No evidence for photon pile- up is found in 
the data. All the spectra are shown in Fig. [3] where their composite nature is clearly visible: 
while at lower energies the thermal spectrum prevails, with the flux values in the 4 different 
extraction regions proportional to EEF, at higher energies almost all the photons detected 
are produced by the PWN and the flux values are proportional to the spatial distribution of 
the nebula surface brightness. Since background estimate will be provided by the spectral 
fitting, we avoid background subtraction. We chose the energy range 0.2 — 10 keV and we 
rebinned the spectra in order to have at least 30 counts per spectral bin and no more than 
3 spectral bins per energy resolution interval. Owing to the extended nature of the diffuse 
nebular emission we generated effective area files with the prescription for extended sources, 
without modelling the PSF distribution of the source counts. 



3.2. EEF From XMM-Newton telescopes PSF 



The XMM-Newton telescopes radially a yeraged PSF profile a t distance r from the aim- 
point is well described by the King's profile (IGondoin et al.lll998l ): 



PSF(r) = A 



1 + 



where the core radius r c and the slope a vary as a function of the photon energy and off-axis 
angle and A is a constant normalisation factor. The calibration files provide the values of r c 
and a, for different photon energies and off- axis angles. 



Since the Vela pulsar thermal spectrum peaks at about 0.5 keV, and the source is 
always at the centre of the field-of-view, we used the parameters r c = 1.59 and a = 5.504, 
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corresponding to an energy of 475 eV and an offset angle of 0.0 degrees. 

The EEF, i.e. the fraction of energy collected within a certain radius R from the source 
distribution centroid, is simply the integral of (jTJ times the radius over the radius, normalised 
at 5': 

i i 

1.59-1 



EEF(R) = J . (2) 



, 5.504' , 



1.59-1 



Since the source is not strictly monochromatic, we calculated the fraction of energy 
encircled in the spectra extraction regions also for 250 eV and 800 eV, in order to obtain a 
"confidence interval" . 



3.3. Surface Brightness Map of the PWN Prom Chandra Data 

We used Chandra data to produce a surface brightness map of the nebular emission. 
In the Chandra /HRC imag es more than 97% of the pulsar counts falls within 3" from the 



centroid of the distribution (IWeisskopf et al.ll2002l . and references therein). Thus, outside this 



small region, the images represent a good approximation of the brightness spatial distribution 
in the energy range 0.1 — 10 keV . 

Since we needed an estimate of the PWN flux in the proximity of the position of the 
pulsar, a two dimensional model of a point-like source flux distribution was obtained with a 
simulation performed with ChaRT^I and MARX0. A 2D fitting was then performed on the 
HRC images with the addition of a spatially constant contribution from the PWN. Best fit 
count rates found for the region between 2.64" and 3.46", corresponding to 20 — 26 pixel, 
are reported in table [2] for the three different HRC-I exposures. Following the criterion 



by iHelfand et al.l (120011 ). the region inside 3.46", 26 pixel, was excised from all images and 



replaced with a poissonian distribution with a mean value equal to the best fit count rate. 

In order to take into account the angular resolution of the XMM telescopes, the maps 
obtained were convolved with a Lorentzian kernel of T — 7.125" corresponding to a FWHM 
of the distribution of 6.6", the nominal FWHM of the XMM mirrors. 

Finally, the convolved image was used to compute the encircled PWN fraction in each 
extraction region of our XMM image. 



6 http : / / cxc . harvard . edu/ chart/ 
7 http : //space . mit . edu/CXC/MARX/ 
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3.4. Background Estimate From Spectral Fitting 

The XMM spectra extracted from the four annuli (Fig. [3]) were fitted simultaneously 
using the following combination: 

(interstellar absorption) x ( (PSR model) + ej(PWN model) ) 

where pi and q are the PSR and PWN encircled fractions within the i th extraction 
region (see Sect. 13.21 and 13.31) . The interstellar absorption does not depend on i. The PWN 
model is a power law whose index 7, is allowed to vary among different regions, while the 
model describing the PSR emission does not depend on i. Our approach simultaneously 
resolve the PSR and PWN emission and compute best fit parameters for both the PWN and 
the PSR models (described in the following sections). 

The spectral fitting was performed with XSPEC v. 11. 3, in the energy range 0.2 — 10 keV. 
The values of the coefficients are given in table [3j and represent the different contribution 
of each region to the total flux of the point-like and diffuse emission. Uncertainties on the 
pulsar as well as on the nebular coefficients arise from the following caveats: 

• the pulsar spectrum is far from monochromatic; 



the XMM-Newton PSF is known within an uncertainty of ~ 2% (IKirschl 120061 ) 



the P WN emission is known to be highly variable over a small time scale (jPavlov et al. 



200l|); 



• the Lorentzian kernel differs from the telescopes PSF; 

» the HRC efficiency differs from the pn one and also the telescopes effective areas have 
different dependencies on the energy of the incoming photons. 

Thus, to account for all uncertainties, an overall systematic error of 5% has been added. 



3.4-1- EPIC/pn Vela Pulsar Wind Nebula spectrum 

The overall (absorption + PSR + PWN) best fit model, summarised in tabled allows us 
to determine the nebular background flux and photon index for all spectra. In particular we 
found that the inner circle contains 7.60 ± 0.04 cts s _1 from the PSR and 2.58 ± 0.02 cts s _1 
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from the PWN (note that only 58% of the pulsar counts are contained in this region - see 
table E] and 0j). 

From table [5] we can see that the non-thermal nebular emission, whose photon index in 
the 4 concentric regions are ji, 72, 73 and 74, become softer as the distance from the pulsar 
increases. The result is co nsistent with Chandra- ACTS spectra obtained from 30 November 



2000 observation (see also iKargaltsev &: Pavlovl 120041 1 . The power-law normalisation yields 



the PWN flux within 40"from the PSR in the EPIC/pn data, its value turns out to be 
Fx = 4.26 ± 0.03 x 10 _11 erg cmT 2 s~ x , in the range 1 — 8 keV. A large part of the 
PWN lies outside such region. In order to estimate the flux of the remaining portion of 
the PWN inside the pn field of view, we extracted a spectrum from the whole detector, 
excluding the 40"circle centered on the PSR, as well as a stripe affected by PSR out-of-time 
events. Ad-hoc response matrix and effective area file were generated. We fitted an absorbed 
power law model to the spectrum in the range 2 — 8 keV, which yielded a best fit photon 
index of 1.61 ± 0.02 and an unabsorbed 1 — 8 keV flux of (1.51 ± 0.02) x 10~ n erg cmT 2 
(^2 _ 0.99, 140 d.o.f.). Thus, the total flux of the PWN inside the pn field of view is 
Fx = (5.77 ± 0.04) x 10 _11 erg cmT 2 s~ x , corresponding to a luminosity of 

L x = (5.74 ± 0.04) x 10 32 erg s" 1 (3) 



at the parallactic distance (~ 287 pc). Such value is in good agreement with IKargaltsev et al 

(SqJ. 



3.5. EPIC/pn Vela pulsar spectrum 



The spectral distribution of the ~ 500 000 Vela PSR photons detected in the inner 
10" radius circle, can be described as the superposition of different components. A simple 
blackbody thermal model does not fit well the observed spectra, that appear harder than a 
simple planckian emission. The structure of the residuals suggests the addition of a second 
"pulsar" component: both a second blackbody or a power-law gave acceptable fits. Similar 
results were obtained with a magnetised hydrogen atmosphere + power-law model. The 
interstellar column density of ~ 2.6 x 10 20 cm" 2 , i nferre d fro m the differ e nt mo dels, is in 
good agreement with the results of iPavlov et al.l (120011 ) and iMori et al.l (12004 ) . For the 



double blackbody model we obtained xl — 11 (1010 degrees of freedom). Assuming a 
parallactic distance of 287 pc, the soft spectrum is described by a cooler (Tbb = 1.06 ± 
0.03 x 10 6 K) component, with a radius i?bb = 5. lio'3 km and a smaller and hotter one 
of T BB = 2.16±g:gf x 10 6 K and .Rbb = 0.73±g;g? km. An equally good fit was obtained 
with a 7 = 3.48+Q Q0 power-law; in this case the blackbody is found to have a temperature 
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T 



1.49lo:°2 x 106 K and a radius R 



2.0^Q4 km. The power-law is much steeper than 



that observed in all other X-ray emitting pulsars. 

A hydrogen atmosphere plus power-law model yields a slightly better fit (xt = 1-0 for 
1010 degrees of freedom). The surface temperature, for a radius R = 10 km and a magnetic 
field B = 10 12 G, is T ha = 0.681 ± 0.004 x 10 6 K; the best fit distance D = 269^ pc agrees 
with the parallattic measure ments. The power-la w photon index 7 = 2.8 ± 0.2 is however 
steeper that the one found by lPavlov et al.l (120011 ). 



In Fig. HI we plotted the V ela unfolded spectrum together with those of Geminga, PSR 
B1055-52 and PSR B0656+14 JPe Luca et al.ll2005h . In the case of Vela we subtracted the 
PWN contribution as estimated from spectral fitting. Given the similarities between the four 
XMM-Newton pulsar spectra, we decided to analyse the Vela data usi ng a phenomenolog- 



ical model encompa ssing two blackbody and a power law component (jCaraveo et al.l 12004 
De Luca et al1l2005h . 



4. Timing Analysis 



The temporal resolution of the pn CCDs in "Small Window" readout mode is 5.6718 
ms. Thus, the time-tagged XMM-pn photons are perfectly suited to study the modulation of 
the Vela PSR flux as a function of its ~ 89 msec rotation period. The photon arrival times 
were corrected for the discrete sampling due to the CCD readout and converted to the solar 
system barycentre with the SAS task barycen. Source counts were extracted from the same 
region used for spectral analysis after subtraction of the non point source contribution. 

Epoch foldi ng yields a strong signal at the pulsar frequency. Following the prescription 
of iLeahyl (119871 ) for high accuracy period determination and error evaluation, we found 
P101 = 0.08933185 ± 0.00000001 s and P 20 i = 0.089331857 ± 0.000000007 s, respectively for 
observation 0111080101 and 0111080201. Such values agree with the contemporary radio 
ephemeris, kindly provided by ATNF (R. Dodson, private communication) and reported in 
table [7J The temporal series were then merged and folded using the radio / and /, to align 
in phase X-ray and radio light curves. 

The overall (0.2-10 keV) X-ray light curve shows 3 broad peaks per period marked 
XSi, XS 2 , XS 3 in Fig. [5J The first and highest peak is phase aligned with the first 7-ray 
peak and follows the radio peak by ~ 0.15 in phase. The second one, the lowest in this 
energy range, reaches its m aximum at cp — 0.45 — 0.50, and co rresponds to the RXTE Peak 
2-Soft Jrlarding et al.l l2002h and UV P2 S jRomani et al.ll2005h . The third peak, which has 
an intermediate intensity, occurs at tp — 0.80 — 0.85; it appears in the soft X-ray light curve 
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and is not present at any other wavelength. The nebular background emission accounts for 
25.5% of the observed flux; the net count rate has a pulsed fraction of 9.2 ± 0.3% (where 
pulsed fraction is defined as the ratio between number of counts above the minimum and 
total number of counts). 

We note the presence of a pulsed signal also at energies above 2 keV, where the contri- 
bution from the two blackbody thermal emission model used to describe Vela spectrum is 
negligible. This would imply a non thermal origin for the pulsed emission above 2 keV, a 
component that was hidden by the nebular emission in the phase-integrated spectrum. We 
estimated the 2-8 keV Vela light curve to be about 100% pulsed. However, we stress that, 
in such an energy range, the source signal accounts for a small fraction (less of 5%) of the 
photons detected, so that a small error in the background estimate would translate into a 
big uncertainty in the pulsed fraction. 

Fig. [6] summarizes the multiwavelength behaviour of the Vela pulsar. Significant dif- 
ferences are seen in the pulse profile emerging from the two energy-resolved EPIC/pn light 
curves: in the energy range 2-8 keV the second peak is stronger than at lower energies and 
the third peak is not observable; the phase interval 0.7 < ip < 1.2 appears more complex at 
energies above 2 keV, with 3 distinct peaks, two coinciding with the optical third and fourth 
peak, and one phase aligned with the first hard X-ray and 7-ray peak. 



5. Phase-resolved Spectral Analysis 

Since in the XMM-Newton domain the pulse shape of the Vela pulsar changes with 
energy, we expect a variation of the spectral shape as a function of the PSR rotation. To study 
such a modulation we extract spectra from 20 different phase intervals. Th e spectra were 



rebinn ed in order to have at least 30 counts per bin. Following the approach of lCaraveo et al. 



koo4 ) and be Luca et all hood ) we have compared the phase-resolved spectra with the two 



blackbody best fit model using the two normalisation coefficients as free parameters. 

Such an approach works well for all the phase resolved spectra but for that encompassing 
phase interval 0.45 < ip < 0.55 which is characterized by high energy residuals, impossible to 
account for with the two blackbody model. The spectrum of such a phase interval is shown 
in Fig. 

Leaving the temperature of the two blackbodies and the hydrogen column density fixed 
to the values found above, a power-law with photon index 2.2^q' 3 is required to fit the 
spectrum (Fig. [8]) for the 0.45 < ip < 0.55 interval. The normalization of the power-law is 
~ 10 times lower than the nebular emission in this phase interval, thus, its contribution to 
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the phase-integrated spectrum is ~ 100 times lower than the nebular one and it cannot be 
seen in the total Vela spectrum. 

Following the detection of a non-thermal component, we decided to perform spectral 
fitting with a three component source model (two blackbodies + power-law) plus the back- 
ground model to account for the nebular emission. As in the previous case, the blackbody 
temperatures where fixed to the values best fitting the total Vela spectrum while a power-law 
photon index of 2.2 was used. The best fitting parameters are given in table [8] while an ani- 



mated version of the phase resolved spectra can be found in http:/ /www. iasf-milano.inaf.it/~deluca/vela/ 



From the blackbody emissions, we can compute blackbody radii (at parallactic distance). 
In Fig. [9] we reported the variation of blackbody radii and power-law normalization as a 
function of the rotational phase. The hotter blackbody presents a low modulation, with a 
single broad peak per period. The cooler one shows a more complex modulation with possibly 
two peaks per period, with a small "dip" between them, well aligned with the radio pulse. 
The cooler component first maximum (R~ 5.3 km) trails the radio pulse phase by about 
10 ms (~ 0.10 — 0.15 in phase). After such maximum, the emitting area suddenly decreases 
to reach its minimum (R~ 4.7 km), at if ~ 0.30, the same phase of the absolute minimum of 
the Vela soft X-ray flux; the cool blackbody radius then slowly grows and reaches its second 
and more pronounced maximum (R~ 5.5 km) at phase 0.8-0.85 responsible for the third 
XMM light curve peak which is not seen in hard X-rays, nor in the ultraviolet or optical 
light curves. 

The hotter blackbody peaks, with a maximum radius of about 0.77 km, at the same 
phase of the first cool blackbody peak, trailing by 0.1-0.15 in phase the radio pulse, long 
known to mark the pulsar polar region. The minimum of the hot blackbody emission occurs 
at about 180° from the maximum and trails the minimum of the cooler emission by ~ 0.3 
in phase. The transition from the minimum to the maximum state is not symmetric, as the 
growth is sharper than the descent. 

The non-thermal emission is present with a narrow peak between if ~ 0.4 and (p ~ 0.6, 
between the two thermal emissions minima. The power-law normalisation is found to be 
consistent with zero in all other phase intervals. Such non thermal component is responsible 
for the second peak observed in the 2-8 keV EPIC light curve and seems to be connected 
with the narrow "spike" in the NUV/FUV light curves and with the "peak 2-Soft" detected 
with RXTE. The pow er-law photon index also agrees with the value found in RXTE phase- 



resolved spectroscopy (IHarding et al.ll2002l ) in the same phase interval. The peak corresponds 



to the "Leading Wing 2" in the EGRET energy range. 
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Discussion 



Considering the different spectral components identified in the EPIC/pn data we note 
that the maximum extension of the hotter thermally emitting region is observed when look- 
ing at the polar region of the neutron star. The radius of this blackbody component 
i s in agreement with the radius of a polar cap within a simple magnetic dipole model 
(IGoldreich fc Julian! Il969l ). For the rotational velo city of the Vela pulsa r , the polar cap 
radius would be R PC = R^RVt/c) 1 ^ 2 = 0.485 km. lHarding fc Muslimovl (120021 ) made an 
estimate of the energetics for an hot spot re-heated by curvature particle downflow: with 
the Vela parameters the bolometric luminosity would be L = 6.01 x 10 31 erg s _1 and the 
temperature T = 2.66 x 10 6 K. Both values agree with what we found. Since the Vela pulsar 
is known to be an inclined rotator, a terrestrial observer would face a single polar region 
during rotation. Assuming that the hot component is produce d by the polar cap region 
heated by return current, the gravitational bending (jPagd Il995l ) would result in a shallow 
modulation of the emitting surface, in agreement with our results. 

The cooler blackbody would represent the radiation from the remaining part of the 
star surface. The radius inferred from the phase-averaged spectrum is too small to fit in 
any proposed equation of state for a star composed mainly of neutron. We also observe a 
~10% modulation of this spectral component as a function of the rotational phase. Mag- 
netospheric reprocessing of the thermal photons emitted from the surface could provide a 
phase-dependent "obscuration" of a fraction of the neutron star surface, depending on mag- 
netic field configuration and viewing geometry. The phenomenon of the mag netospheric 
"blanket", e.g. cyclotron reso nance scattering by pla s ma a t a few stellar radii (IRuderman 



20031 ). originally proposed by lHalpern fc Rudermanl (119931 ) as an explanation of the soft 
thermal emission of Geminga, could provide the physical basis for the observation of a 
phase-dependent emitting area. 

Anisotropic heat transfer from the interior of the neutron star, would also provide a 
surfa ce temperature far from un iform, with the polar regions hotter than the equatorial 
ones (IGreenstein fc Hartkdll983l ). Such non-uniform temperature distribution would result 
into a modulated X-ray thermal flux. Alternatively, using a uniform su rface t empe rature 
approximation, this translates into a modulation of the emitting surfaces. iPagd (119951 ) show 
that in such a case a modulation of a few percent is expected, owing to gravitational light 
bending. However such a modulation would have to be phase aligned with that of the hotter 
component. The complicate modulation we observe for the cool blackbody component is 
not easy to reconcile with such a picture, pointing to a more complex surface temperature 
distribution, and/or to magnetospheric reprocessing as discussed above. 



The non-thermal component detected in a small phase interval 180° from the radio 
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pulse, when the thermal components are at their minimum, might be produced by particles 
accelerated in the star magnetosphere. 



Conclusions 



The spectral analysis of the XMM-Newton data of the Vela pulsar confirms the thermal 
nature of its emission pointing towards an X-ray phenomenology similar to that seen in 
older neutron stars. Both the overall spectral shape and the phase-resolved behaviour are 
reminiscent of what has been found for middle aged pulsars such as Geminga, PSR B0656+14 
and PSR B1055-52, affectio nately called the "Three Musketeers" on view of their similarities 
( jBecker fc Truemperl Il997l ). Indeed, the Three Musketeers could be described, even with 
some caveats, with a simple common phenomenological model in which the differences in 
the variations of the "hot spot" emitting cLFGcLS clS cL function of the rotational phase could be 
ascribed to the different viewing geometries of the three neutron stars which are known to 
be different: orthogonal rotator seen perpendic ular to the rotationa l axis for PSR B1055-52 
and almost aligned for PSR B0656+14 (see e.g. iDe Luca et al.ll2005l for a full description of 
such a picture.). 

While Vela is certainly hotter than the older musketeers, its overall blackbody emitting 
radius is rather small leading to a low luminosity , well below the value expected in the stan- 
dard cooling scenario described by lTsuruta et al.l (120021 ). Using a hydrogen atmosphere eases 
the radius problem but, yielding lower temperature, does not change the source luminosity 
which rests on a precisely determined parallactic distance. 

Moreover, although Vela is known to be an inclined rotator (the angle between the 
magnetic dipole vector and the rota tion axis is a ~ 60° , while the line of sight is inclined 
with respect to the latter by £ ~ 55°, iHelfand et aJj|200l| , and references therein), our phase 
resolved spectroscopy points to a surprisingly small modulation, typical of aligned rotators, 
such as PSR B0656+14, rather than orthogonal ones. 

The presence of a non-thermal component sharply pulsed in antiphase with the radio 
pulse is a surprising result of our analysis and could be exploited to better constrain the 
geometry of the emitting regions with the goal to unveil the emission mechanism behind 
Vela complex multiwavelength behaviour. 

The remarkable body of evidence collected so far on the Vela pulsar has been ge nerally 
interpreted in the framework of the outer gap model (jRomani &: Yadigaroglulll995l ) where 
e + / - couples fill the pulsar magnetosphere high above the star surface. Indeed, such a model 
predicts an X-ray non-thermal broad peak at (p ~ 0.35 — 0.5, a feature clearly seen both in 
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RXTE and XMM light curves. 

Recently, the Vela pulsar pheno menology has been accounted for also by the two pole 



caustic model (IDyks fe Rudakll2003l ) to supply the high energy particles. 



XMM-Newton findings provide a useful link between the star thermally emitting surface 
and its particle filled magnetosphere. Such two-regime, phase-resolved view probes the star 
geometry as it rotates and provides inputs for future tri-dimensional realistic models of the 
pulsar magnetosphere. 

We thank R.Dodson for providing Vela pulsar radio ephemeris simultaneous with the 
XMM-Newton observations. We thank A.Possenti for his help with the radio data. We 
would also thank A.Harding and M.Strickman for their help to understand the RXTE phase 
alignment. XMM-Newton and Chandra data analysis is supported by the Italian Space 
Agency (ASI) under contract ASI/INAF 1/023/05/0. ADL acknowledges an ASI fellowship. 
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Table 1: Journal of XMM- Newton observations. 



Obs ID 


Date 


MJD 


Obs. Time (ks) 


Instrument (mode) 


Good Time (s) 


Live Time (s) 


0111080101 


2000 Dec 1 


51879 


41.11 


pn(SW) 


37 590 


26 366 










MOSl(SW) 


37648 


37 241 










MOS2(LW) 


37694 


36 597 


0111080201 


2000 Dec 2 


51880 


61.81 


pn(SW) 


50 887 


35 693 










MOSl(SW) 


49 067 


48 524 










MOS2(LW) 


50410 


48 940 



Note. — For the two different XMM datasets are listed the start date of the observation in standard 
UTC time and in Modified Julian Day (MJD), the total observation time, the instrument that performed 
the observation and the operational mode (SW= "Small Window", LW= "Large Window"), the good time 
intervals after removing the Soft Proton Flares, the Live Time after correction for charge transferring and 
readout time. 



Table 2: 2D best fit counts in the annular region 2.64"-3.43"for the three HRC-I images. 



Obs. Id 


364 


1518 


1966 


cts pixel -2 


1.43 ±0.04 


1.41 ±0.04 


1.23 ±0.04 


cts pixel -2 s -1 


0.0017 


0.0016 


0.0014 



Table 3: Total number of counts. 



Extraction 


Total number 


PSR EEF 


PWN flux 


region 


of photons 


(%) 


(counts) 


10" 


671 021 


0.582342 


17729.6985 


10"-20" 


541 089 


0.215818 


40962.9075 


20"-30" 


431813 


0.0770648 


42880.3005 


30"-40" 


293 991 


0.0377001 


29622.9375 



Note. — Total number of photons collected by the EPIC/pn in the different spectra extraction regions. 
Also shown are XMM- Newton EEF and the PWN counts estimated on the Chandra/HRC surface brightness 
map (degraded to the XMM angular resolution). The PWN parameters e» used in the fit are proportional 
to such counts. We assumed ei=l. These values have been used to compute the PSR and PWN coefficients 
for the spectral fitting. 
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Table 4: Pulsar and Nebula contributions. 



Extraction 


PSR flux 


PWN flux 


EPIC/pn counts 


region 


(counts) 


(counts) 


(0.2-10 keV) 


10" 


471592 


160114 


631711 


10"-20" 


164272 


355 540 


519824 


20"-30" 


54 544 


363174 


417689 


30"-40" 


36491 


247743 


284238 



Note. — Number of photons collected by the EPIC/pn in the different spectra extraction regions, emitted 
from the Vela pulsar and from its PWN, as derived from best fit spectral model fluxes. 
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Table 5: Best fit parameter for double blackbody model. 



Parameter Best fit value 

Cool blackbody (bb) + Hot blackbody (BB) 
nH ( x 10 22 cm- 2 ) (2.59 ± 0.01) x 1(T 2 
k B T hh (keV) . 0.091 ± 0.003 

norm hh 18000^^ 

k B T BB (keV) 0.186±g;g°i 
norm BB .... 380±?g 

71 J-^O' -0.007 

ry n i 4qc;+o.oo6 

72 i -^ Jd -0.006 

73 i - OJJ -0.006 

-0.007 
-0.006 



74 1.651+ 0007 

normp WN . . . (1.033 ± 0.006) x 10 



-3 



Note. — Best fit parameter for the PSR double blackbody model and for the PWN. Spectra from the four 
regions described in Sect. l3.1l have been used. nH is computed wit h a photo-electric absorption model (model 
wa bs in XSPEC) using Wisco nsin ( Morrison fc McCammonlll983h cross-sections and the relative abundances 
bv lAnders &: Ebiharal (j 19821 ) . Blackbody normalizations are defined as (R/Dio) 2 , where R is the emitting 
radius in km and Dio is the distance in units of 10 kpc. The PWN photon indices 71, 72, 73 and 74 are the 
best fit values for the four extraction regions, from the innermost circle to the outermost annulus (see text). 
As explained in the text, a common PWN normalization parameter is simultaneously fitted to all spectra, 
since the PWN flux fraction in different region is computed a priori using the Chandra surface brightness 
map (degraded at the XMM angular resolution) and is kept fixed during the fit. The normalization of the 
PWN in the i th region is obtained by multiplying the value quoted in the table (in units of photons cm" 2 s~ x 
keV -1 at 1 keV) by the PWN parameter (see text and table [3]). Errors are computed at a 90% confidence 
level for a single interesting parameter. 
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Table 6: Physical characteristics of the vela pulsar 



Parameter Best fit value 

Cool blackbody (bb) + Hot blackbody (BB) 



T bb (K) ... (1.06 ±0.03) x 10 6 

R hh (km) . 5.06j£g 
L bb (erg s- 1 ) (2.3+°i) x 10 32 



T BB (K)... (2.161°,$) x 10 

R BB (km) . 0.73l°;°? 
L BB (erg s" 1 ) (8.3l 2 ; 2 ) x 10 31 



L toi (erg s -i) (3.11° j) x 10 32 

l pwn ^ erg g -i)i ( 5 74 ± 0>04 ) x 1Q 32 



Note. — Physical characteristics of the Vela pulsar as inferred from the cool blackbody (bb) plus hot 
blackbody (BB) best fit model (see table We assumed 58% of the energy emitted encircled in the 10" 
radius region. Emitting radii are computed for a distance of 287 pc. The bolometric blackbody luminosities 
(-Lbb) Lbb and L^ 1 ) are computed from Stefan-Boltzmann law, for best fit temperature and radius. The 
PWN luminosity (Lx WN ) is computed for the energy range 1-8 keV at the parallactic distance, using the 
flux measured within the entire pn field of view (see Sect l3.4.Tj) . 



Table 7: Vela radio ephemeris (R. Dodson, private communication). 



t in MJD . . . 

/ (Hz) 

p (s) 

/(Hzs- 1 ).... 
P (s s" 1 ) 



51881.000000049 
11.1942146219182 
0.0893318588016 
-1.56297 x 10~ n 
1.24728 x 10~ 13 



Table 8: Results of phase resolved spectroscopy 



Phase 


bb radius 


BB radius PL norm xl (d.o.f.) 




(km) 


(km) (7 s -1 cm -2 keV -1 ) 






00- 





05 


5 


12 


± 





06 


0.73 ± 0.02 


(1.6 ± 1.2) x 10~ 4 


1 


00 


(136) 





05- 





10 


5 


22 


± 





06 


0.75 ± 0.01 








99 


(137) 





10- 





15 


5 


33 


± 





06 


u - ' ' -0.02 





1 


10 


(137) 





15- 





20 


5 


27 


± 





06 


0-77+°;° 2 


(1.9 ± 1.2) x 10~ 4 


1 


10 


(143) 





20- 





25 


5 


11 


± 





06 


0.75 ± 0.01 





1 


28 


(137) 





25- 





30 


4 


80 


± 





06 


0.76 ± 0.01 





1 


41 


(133) 





30- 





35 


4 


TO 

73 


1 

± 





06 


°- 74 -0.02 





1 


10 


/ 1 O A \ 

(134) 





35- 





40 


4 


84 


± 





06 


7S+ - 01 

u - ' °-0.02 





1 


51 


(138) 





40- 





45 


4 


88 


± 





06 


0.73 ±0.02 


(1.8 ± 1.2) x 10" 4 


1 


17 


(138) 





45- 





50 


4 


88 


± 





06 


0.72 ±0.02 


(5.7 ± 1.2) x 10~ 4 





84 


(140) 





50- 





55 


4 


87 


± 





06 


0.69 ±0.02 


(5.0 ± 1.2) x 10~ 4 


1 


13 


(144) 





55- 





60 


5 


08 


± 





06 


66 +0 - 02 


(2.0 ± 1.2) x 10~ 4 


1 


03 


(139) 





60- 





65 


5 


09 


± 





06 


6 c ; +0 - 01 

u - oo -0.02 





1 


08 


(138) 





65- 





70 


5 


13 


± 





06 


0.71 ±0.01 





1 


08 


(137) 





70- 





75 


5 


18 


± 





06 


0.75 ±0.01 





1 


26 


(136) 





75- 





80 


5 


38 


± 





06 


0.73 ±0.01 





1 


08 


(134) 





80- 





85 


5 


46 


± 





06 


0.74 ±0.01 





1 


09 


(133) 





85- 





90 


5 


33 


± 





06 


0.71 ±0.02 


(5.0+1°) x 10 - 5 


1 


27 


(134) 





90- 





95 


5 


16 


± 





06 


72+ - 01 

u - ' z -0.02 


(0.3+ ° 3 ) x 10- 5 


1 


14 


(136) 





95- 


1 


00 


5 


05 


± 





06 


0.75 ±0.01 





1 


11 


(134) 



Note. — Cool blackbody (bb), hot blackbody (BB) radius and power-law normalisation as function of 
the rotational phase. Temperatures of both blackbody components were fixed at the best fit values for the 
phase-integrated spectrum. The power law photon index was left as a free parameter in the 0.45-0.55 phase 
interval (best fit value of 2.27_q; 4 ). For the remaining part of the cycle the photon index was fixed at the 
best fit value for the 0.45-0.55 phase interval. For the radius computation we assumed 58% of the energy 
emitted encircled in the 10"radius region and a parallattic distance of 287 pc. Errors are computed at a 90% 
confidence level for a single interesting parameter. 
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Fig. 1. — EPIC/pn image of the Vela pulsar. The cross marks the radio position and the 
box marks the region covered by Chandra image of Fig. [2J The spectral extraction regions 
(Sect. 13. ip are also marked. 
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Fig. 2. — Chandra-HKC image of the Vela pulsar. 
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Fig. 3. — Observed EPIC/pn spectra from the 4 different extraction regions. 
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Fig. 4. — XMM-Newton spectra of the Three Musketeers and the Vela pulsar (IDe Luca et al. 
20051 ) . The Vela pulsar thermal spectrum drops below the harder nebular emission at energies 
above ~ 2 keV. 




0.5 1 1.5 2 

Phase 



Fig. 5. — Folded light curve of the Vela pulsar in the 0.2-10.0 keV energy range as observed 
with the EPIC/pn camera. Two rotational periods are represented for clarity. Phase 
corresponds to the radio peak. 
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Fig. 6. — Multiwavel ength light curves of the Vela pulsar from radio to 7-rays (lHarding et al. 
20021 ; lKanbach!l2002l ). Two rotational periods are represented for clarity. In the EPIC/pn 
light curves, error bars are also plotted. 
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Fig. 7. — Observed spectrum for the phase interval 0.45-0.55 and best fit double blackbody 
model folded with the instrument response. The two blackbody components of the pulsar 
spectrum are shown (green curve: cool blackbody; red curve: hot blackbody) together with 
the nebular contribution (yellow curve). The lower panel shows the residuals in unit of 
standard deviation. Above ~ 2 keV, where the nebular contribution is overwhelming, the 
data are systematically above the best fit encompassing the double blackbody together with 
the nebular contribution. 
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Fig. 8. — As for Fig. [7] but with a double blackbody plus power-law model for the pulsar 
emission. Note the reduction of the high energy residuals. 
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Fig. 9. — Variation of the three spectral component (cool and hot blackbody radius and 
power-law normalization) as a function of the rotational phase. Two rotational periods are 
plotted for clarity. 



